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Abstract. We consider the question of cross-correlation measurements using 
Virgo and the LSC Interferometers (LIGO Livingston, LIGO Hanford, and 
GEO600) to search for a stochastic gravitational-wave background. We find that 
inclusion of Virgo into the network will substantially improve the sensitivity to 
correlations above 200 Hz if all detectors are operating at their design sensitivity. 
This is illustrated using a simulated isotropic stochastic background signal, 
generated with an astrophysically-motivated spectral shape, injected into 24 hours 
of simulated noise for the LIGO and Virgo interferometers. 
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1. Introduction 

There are four kilometre-scale interferometric gravitational-wave (GW) detectors 
currently in operation: the 4 km and 2 km interferometers (IFOs) at the LIGO Hanford 
Observatory (LHO) (known respectively as HI and H2), the 4km IFO at the LIGO 
Livingston Observatory (LLO) (known as LI), and the 3km Virgo IFO (known as 
VI). The LLO and LHO IFOs are currently conducting operations at their design 
sensitivity, and operate together with the 600 m GEO600 IFO (known as Gl) under 
the auspices of the LIGO Scientific Collaboration (LSC). Virgo's first full science run 
commences in May 2007. 

One of the signals targeted by ground-based GW IFOs is a stochastic GW 
background (SGWB), which can be either of cosmological or astrophysical origin, in 
the latter case being produced by a superposition of unresolved sources. The standard 
technique to search for a SGWB looks for correlations in the outputs of multiple 
detectors. We describe in this paper how the inclusion of correlation measurements 
involving Virgo could improve the sensitivity of the current LLO-LHO network. 
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HL HV LV GH GL GV 



712(0) -0.89 -0.02 -0.25 0.42 -0.32 -0.08 

(27TT12)- 1 (Hz) 15.9 5.8 6.0 6.3 6.3 49.8 
7^ v (Hz) -1.7 5.0 -5.9 3.0 4.3 -17.2 

Table 1. Limiting behaviour of overlap reduction functions of detector pairs. H 
refers to either of the IFOs at the LHO site, L to LLO, V to Virgo, and G to 
GEO600. At / = 0, the ORF is determined by the alignment of the detectors. 
The reduced inverse light travel time (2ttTi2)~ 1 gives a characteristic frequency 
for the onset of "high-frequency" behaviour, which includes a sine function of the 
ratio of / to that frequency. However, the limiting form in {2) also includes a 
geometric projection factor, leading to an overall envelope 7™ v which is shown in 
the third row. In particular, while the light-travel time Thl is l ess than Tuv or 
T^v, the projection factor more than makes up for this, which makes the mean 
amplitudes of 7hv(/) an d 7Lv(/) a * high frequencies larger than that for 7jjl(/)- 



2. All-sky sensitivity at design 

2.1. Observing Geometry 

The effect of a SGWB is to generate correlations in the outputs of a pair of GW 
detectors, which can be described for an isotropic background in the Fourier domain 

by 

(K(f)h 2 (f)) = ^(/-/')7i 2 (/)5 gw (/) . (1) 

The raw correlation depends on the (one-sided) power spectral density S l gw (/) the 
SGWB would generate in an IFO with perpendicular arms, as well as the observing 
geometry. The geometrical dependence manifests itself via the overlap reduction 
function (ORF)[T], which can be written as[2] 

112(f) = d lab df i- JJ d 2 n A e «*/Mft-ft)/c (2 ) 

where each IFO's geometry is described by a response tensor constructed from unit 
vectors x and y down the two arms 

d ab = l -(x a x b -y a V b ) , (3) 

fi,2 is the respective interferometer's location and P TTri cd is a projector onto traceless 
symmetric tensors transverse to the unit vector n. At zero frequency, the ORF is 
determined entirely by detector orientation. The LHO and LLO sites are aligned as 
nearly as possible given their separation on the globe, so that Jhl{0) — —0.89. In 
contrast, the Virgo and GEO600 sites are poorly oriented with respect to one another, 
so 7gv(0) = —0.08. However, the frequency-dependence of the ORFs means that the 
situation is quite different at frequencies above 40 Hz, where the IFOs are sensitive. 
In particular, the amplitude of jGv(f) does not drop appreciably for / below about 
350 Hz. For the other pairs, the behaviour is determined by the high-frequency limiting 
form of the ORF, 

712 (/) — 5di Qb P TTSl2a c b d df sinc(2^/T 12 ) = M- sin(27r/T 12 ) (4) 

where T12 is the light travel time between the detector sites and §12 is a unit vector 
pointing from one site to the other. While trans- Atlantic light travel times like Thv 
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Figure 1. The overlap reduction functions for pairs of detector sites. Note 
that the ORF for the two LIGO sites goes off the scale of this plot at 40 Hz, 
which is the "seismic wall" below which LIGO data are too noisy to be of 
any use. The proximity and alignment of LLO and LHO makes HL the most 
favourable pair of detector sites for observations below 150 Hz or so. However, the 
proximity of the Virgo and GEO600 sites means the GV ORF is substantial out 
to higher frequencies, overcoming the low- frequency suppression due to their poor 
alignment. On the other hand, as shown in Table [T]the HV and LV ORFs, while 
they oscillate rapidly with increasing frequency, do not decay as precipitously 
as the HL ORF, making them more favourable than HL (but less than GV) for 
/ > 200 Hz. 



and Tlv are greater than Thl, leading to more oscillations in the ORF, the overall 
envelope 7®f v = {5di a bP TTsi2 cd^2 d )/(^ 7T ^'2) includes geometric projection factors, 
which more than make up for this discrepancy, as summarised in Table [TJ The result 
is that in the full overlap reduction function (Fig. [T]) the typical amplitudes |7hv(/)| 
and |tlv(/)| are larger than the typical |thl(/)| for / > 200 Hz. 



2. 2. Definition of Sensitivity 

The strength of an isotropic stochastic background can be written in terms of the one- 
sided power spectral density S gw (f) it would generate in an IFO with perpendicular 
arms. 

The standard cross-correlation method seeks to measure the amplitude Sr = 
5gw(/)/5(/) of a background whose S gw (f) is assumed to have a specified shape 
<S(/). Given coincident data between times t\ and t<i from detectors with one-sided 
noise power spectral densities (PSDs) P%,2(f), we can make an optimally-filtered cross- 
correlation statistic 

Y= [ [ Sl (t)Q(t-t')s 2 (t)dtdt' (5) 
Jti Jti 

with the optimal filter defined by its Fourier transform 

Q(f)=M 112if)S{f) (6) 
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Figure 2. Sensitivity integrands 1(f) for pairs of detector sites, where pairs 
like HL include e.g., the combined sensitivity of Hl-Ll and H2-L1. To get an 
overall sensitivity, these need to be multiplied by the square of the shape of 
the SGWB spectrum Sg W (/)- The HL pair is the most sensitive for fiat and 
low-frequency spectra, but as shown in the closeup at the right, above 200 Hz, 
pairs involving Virgo are more sensitive when all detectors are operating at their 
design sensitivities. The GH and GL pairs are not shown, since the GEO600 noise 
spectrum [5] means that only the GV pair contributes significantly to the overall 
sensitivity at these frequencies. 



and M chosen so that (Y) = Sr. If the geometric mean of the noise PSDs is large 
compared to S gvr (f), the expected variance of the statistic will be 

-l 




[7i2(/)S(/)] 2 



(7) 



where T = t% — t\ is the duration of the data analysed. The squared signal-to-noise 
ratio of the standard cross-correlation statistic will thus be 

ly\2 poo 

SNR 2 :=\—L = 2TS 2 R [S(f)] 2 l l2 (f) df , (8) 
° Jo 

where we have defined a "sensitivity integrand" which illustrates the contribution to 
the sensitivity of different frequencies: 

Il2(/) W)m (9) 

We plot 1(f) for several pairs of detectors in Figure [2j using the design sensitivities 
in[3l 01 [5]. As shown in [6], the optimal method for combining correlation 
measurements from different detector pairs is the same as that for combining 
measurements from different times: average the point estimates Y with a relative 
weighting of rr~ 2 , and the resulting variance will be the inverse of the sum of the a~ 2 
values. This produces a sensitivity integrand which is the sum of the integrands for 
individual pairs: 

X(/)=^X pair (/) . (10) 

pair 

An immediate application of this is to define sensitivity integrands that combine pairs 
involving HI and H2, e.g., 2"hl = 2"hi,li +2"h2,li- This is the same as using the 
spectrum of an optimally combined H pseudo-detector as described in [7). 
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Figure 3. Combined sensitivity integrands for networks of detectors. In each 
case a network including "H" includes correlations between both HI and H2 and 
the detectors at other sites. As the closeup on the right shows, addition of the LV, 
HV and GV pairs to the HL network increases sensitivity to backgrounds with 
significant power above 200 Hz. 



Table 2. Smallest detectable band-limited background using each of the detector 
networks defined in Fig. [3] In each case, this is the strain power spectrum, in units 
of 10 — 48 Hz -1 , that could be detected with 5% false alarm and 5% false dismissal 
rates, using one year of coincident data at design sensitivity. 

Band H-L H-L-G H-L-V H-L-V-G 

200-300 Hz 5.79 5.43 3.44 3.04 
300-400Hz 18.57 15.37 7.92 5.88 

Figure [3] shows the combined sensitivity for four networks of detectors operating 
at design sensitivity: the existing H-L network, an H-L-G network in which GEO600 
is also operating at design sensitivity, and H-L-V and H-L-V-G networks which also 
include a design-sensitivity Virgo. The H1-H2 pair is not included in these networks, 
because the presence of correlated environmental noise necessitates special treatment 
of this pair [8]. 

Since the power in the faintest detectable "white" stochastic background is 
proportional to the square root of the area under the sensitivity integrand, we see that 
the addition of Virgo to the LHO-LLO network would be most useful in improving 
sensitivity to a narrow-band background peaked above 200 Hz, or to one whose 
spectrum rises with increasing frequency. As an illustration, Table [2] shows for several 
detector networks the faintest detectable background with constant 5 gw (/) in a 100 Hz 
band, assuming one year of observation time and an SNR threshold of 3.29, associated 
with 5% false alarm and false dismissal rates. 

3. Simulations 

To test cross-correlation analyses of LIGO and Virgo data, we injected a simulated 
SGWB signal into simulated LIGO and Virgo noise. The simulated noise data were 
the 24 hours of HI, H2, LI, and VI data, all at nominal design sensitivity, initially 
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Figure 4. The magnetar spectrum used to generate simulated signals, and 
the associated integrand for squared signal-to-noise-ratio. Note that since 
Sgw(/) increases approximately linearly with frequency up to about 400 Hz, this 
spectrum, while broad-band, tends to favour the higher frequencies where LIGO- 
Virgo detector pairs are more sensitive. Beyond 400 Hz, the spectrum is no longer 
linear, but we show the SNR that would result from attempting to detect it with 
a S(f) oc / filter; the corrections are negligible below about 500 Hz, and still small 
throughout the frequency range displayed. Note that integrating the area under 
the H-L-V curve on the right still gives an SNR below 10" 6 from 24 hours of 
data, so we scale up the injected strain signal by a factor of several thousand in 
the simulations described in Section l4l 

generated for simulated searches for GW bursts and inspiralling compact object 
binaries [9,10, known as "project lb" . We chose a spectral shape designed to highlight 
the performance of the LV and HV pairs at / > 200 Hz, but which corresponds to 
a model of an astrophysical SGWB. The spectrum we used is associated with the 
superposition of the tri-axial emission from the extra-galactic population of spinning 
magnetars with type I superconducting interior, as described in model B of [11) . but 
updated using the star formation history of [12] . 

In Figure [4] we show the spectrum and the associated sensitivity integrand in the 
corresponding detectors. Since this spectrum rises with increasing frequency up to 
about 400 Hz, it is useful for illustrating the utility of a network involving Virgo to 
search for a broad-band astrophysical source. Since the model signal would be too 
weak to detect with first-generation interferometers, we scale up the signal strength, 
injecting a signal with the same spectral shape, but a much larger amplitude, in our 
simulations. 

3.1. Simulation algorithm 

The problem of simulation of the signal in a pair of detectors due to an isotropic and 
Gaussian SGWB has been considered previouslyin e.g., [6l[T3l[14]. For this work we 
generalise that to a network of N GW detectors LB We need to satisfy (Q]) for each pair 
of detectors; treating as the elements of a column vector h(/) and {7as(/)} 

I In our case, to simulate signals in HI, H2, LI, and VI, N = 3, because HI and H2 have the same 
response tensor and therefore the same simulated GW signal can be used for both of them. 
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as the elements of a real, symmetric matrisd t(/)j we can write this as a matrix 
equation 

<h(/)h(/') t ) = ^gw(/)7(/W-/') (11) 
If we can define a matrix (3(f) which factors f(f): 

7 (/)=/3(/)/3(/)t, (12) 
then we can generate N independent white noise data streams {?7a(/)} which satisfy 

m)v(n f ) = is(f~n (is) 

and then convert them into the desired coloured correlated data streams via 



Sgw(f) 



h(/) = ^W)r/(/) (14) 
For a given *y(f), there are different choices of rj(f) which achieve the factorisation 

G3- 

Since (|11[) is a covariance matrix, it is positive semi-definite, from which it follows 
(since 5 gw (/) > 0) that j(f) is positive semidefinite as well. -y(f) could have one or 
more zero eigenvalues in the presence of linear dependence between detector outputs. 
A practical example of this is two detectors sharing the same geometry and location, 
such as HI and H2, are included in the network|Jj] We avoid this problem by generating 
a simulated signal for HI and injecting it into both HI and H2. 

In the generic case where ~f(f) is positive definite, we can make the straightforward 
choice of the Cholesky decompositionfTB] . in which [3(f) is a lower diagonal matrix. 
In the case N — 2, this reduces to the form used in e.g. ,[B]. For N = 3, if the diagonal 
elements of ~f(f) are unitjfjj, the explicit form is 

10 \ 

(15) 



W) 



V / l 3 7f2 

723-712713 / 1 + 2712713723 -7l2-7l 3 -723 



However in practice we can simply use a fast iterative algorithm for the Cholesky 
decomposition. 

Other factorisation strategies which treat the different detectors more 
symmetrically (e.g., defining [3(f ) = A(f) 1 / 2 X5(f) where A(/) is the diagonal matrix 
of eigenvalues of ~f(f) and U(/) is the matrix constructed from the corresponding 
eigenvectors) may be more demanding in terms of computational power, but more 
numerically stable when correlations between the detectors are large and off-diagonal 
elements of 'j(f ) are comparable to unity. In particular, this strategy can deal directly 
with the case when f(f) has one or more zero eigenvalues. 

§ For non-isotropic backgrounds, the ORF is complex rather than real, and more care must be taken 
with the definition of the Hermitian matrix f(f). 

|| Less trivial examples can be constructed, for example three detectors in the same location and in 
the same plane. 

This is the case for interferometers with perpendicular arms, but not for GEO600 or for resonant 
bar detectors; see [2] 
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Figure 5. The overlap-and-add algorithm 
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Filtering strategy 

The continuous frequency-domain idealisation (|14|) needs to be applied with care to 
finite stretches of real detector data. In the time domain, the multiplication (fT4| 
amounts to a convolution 

/oo 
K{t-t')r]{t')dt' (16) 
-OO 

with a kernel which is the inverse Fourier transform of 

K(/) = ^43(/) . (17) 

If the time-domain kernel K(r) is negligible outside the interval — t_ < t < t + , we 
can use the standard overlap-and-add strategy to generate a continuous stream of 
time-series data, as illustrated in Figure 

(i) Generate a sequence of "buffers" of white noise data for each of the N detectors, 
each of length T. 

(ii) Convolve each buffer with the kernel K(r) to obtain a time series of length 
r_ + T + t + , with an associated start time r_ before the start and end time t + 
after the end of the original buffer. (This is most naturally done in the frequency 
domain, zero-padding the white noise by r_ at the beginning and r + at the 
end, then Fourier transforming and multiplying by K(/) before inverse-Fourier- 
transforming.) 

(iii) Add together the processed data buffers, overlapping by r+ on one end and r_ 
on the next, producing correlated coloured time-series data of duration T times 
the number of buffers, plus transients of r_ at the beginning and r+ at the end, 
which are discarded. 

This strategy was implemented in code based on Virgo's "Noise Analysis Package" 
(nap) [T5]. 

4. Analysis of simulated data 

The continuous signals described in Section [3] were injected into the "project lb" 
simulated noise and the resulting time series analysed using the MATAPPS stochastic 
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Figure 6. The individual point estimates and one-sigma error bars with signals 
injected at two different levels. The strength of the injection is shown on the 
vertical axis, as a strain power at 200 Hz, while the estimate of this quantity from 
each detector pair is on the horizontal axis. The error bars are drawn vertically 
for ease of reading, and because they are frequentist error bars as quoted. Note 
that the signals are seen at comparable strength in all five pairs of detectors. 



analysis code developed by the LSC. [17] • In particular, the cross-correlation ([5]) 
was performed in the frequency domain without any need to resample the time-series 
data, using a variation of the method described in [18j and applied in |14j . Several 
simulation runs were performed in which same set of simulated signals were injected 
into the four data streams, scaled up by a different factor for each run. The results 
of two of those simulations are shown here. In Figure [5] we plot the individual point 
estimates and error bars in each of the five detector pairs (every combination except 
for H1-H2). The amplitude measure Sr quoted is the one-sided strain PSD at 200 Hz. 
The optimal filter used assumed the shape <S(/) oc /, although the injected signals 
only had that form below about 400 Hz. The analysis was done over a frequency band 
50-500 Hz to avoid difficulties arising from this mismatch, as illustrated in Figured) 
The optimal combination of these results is shown in Table [3l both for the full 
network and for the network consisting only of the two LHO-LLO pairs. Our error bars 
are reduced by 15% via the inclusion of the LIGO- Virgo pairs. Note that although this 
source has more power at the intermediate frequencies favoured by the LIGO- Virgo 
pairs, it is still broad-band. A narrow-banded source with most of its power above 
200 Hz would favour LIGO- Virgo pairs to an even greater extent. 

5. Conclusions and Outlook 

We have demonstrated how the inclusion of LIGO- Virgo and possibly GEO600- Virgo 
detector pairs can enhance the sensitivity of the global GW detector network to 
an isotropic background of gravitational waves, particularly at frequencies above 
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S gw (200Hz) (lO- 48 ^ 1 ) 
Injected H-L Result H-L-V Result 



1.35 
2.70 
6.74 
13.49 



0.39 ±0.98 
1.69 ±0.99 
5.63 ± 1.02 
12.22 ± 1.06 



0.43 ±0.82 
2.31 ±0.67 
6.29 ±0.69 
12.35 ±0.88 



Table 3. The values of .Sgw(200Hz) calculated from the simulated project lb 
data, with associated one-sigma error bars, for the H-L network consisting of the 
Hl-Ll and H2-L1 pairs, and for the H-L-V network consisting of those plus the 
Hl-Vl, H2-V1, and Ll-Vl pairs. We see that including Virgo in the network 
reduces our error bars for this broad-band astrophysical spectrum by 15%. 



200 Hz. As a practical illustration, we have adapted and applied pipelines for 
generating correlated simulated signals in the LSC and Virgo detectors, and for 
analysing coincident data via the standard cross-correlation technique. The specific 
astrophysical model we used (which was chosen because its frequency spectrum was 
peaked at frequencies where LIGO- Virgo pairs at design will be more sensitive than 
LLO-LHO pairs) had to have its amplitude increased to be detectable by any pair of 
first-generation IFOs. Nonetheless, the exercise illustrates how multiple detector pairs 
can be used to discover an "unexpected" background. 

Virgo is not yet at its nominal design sensitivity, but has improved its sensitivity 
markedly over the past year, and its first full science run starts in May 2007, to be 
analysed in conjunction with the end of LIGO's S5 run. 
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